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ABSTRACT 



The simplest partition function, associated with homogeneous symmetric forms S of degree r in n variables, 
is integral discriminant J n \ r {S) = J e~ s ( xl <"'< Xn > dx\ . . . dx n - Actually, 5-dependence remains the same if e~ s 
in the integrand is substituted by arbitrary function f(S), i.e. integral discriminant is a characteristic of the 
form S itself, and not of the averaging procedure. The aim of the present paper is to calculate J n \ r in a number 
of non-Gaussian cases. Using Ward identities - linear differential equations, satisfied by integral discriminants 
- we calculate J213, J214, >^2|5 an d J313. In all these examples, integral discriminant appears to be a generalized 
hypergeometric function. It depends on several SL(n) invariants of S, with essential singularities controlled 
by the ordinary algebraic discriminant of S. 



Contents 



1 Introduction 1 

2 Ward identities 3 

3 Integral Discriminants 5 

3.1 The case of J 2 | 3 12 

3.2 The case of J 2 |4 13 

3.3 The case of J 2 \5 18 

3.4 The case of J 3 | 3 28 

4 Conclusion 33 



1 Introduction 

Averaging with exponential weight 

(0) = / <j)(x)e- s{x) dx 



is an important operation in statistical and quantum physics. Function S(x), which determines the weight, is 
called action. The integration domain of x- variables is a linear space, which is usually infinite dimensional in 
real physical applications: for example, a space of paths in quantum mechanics or a space of field configurations 
in quantum field theory. Infinite dimension of the space introduces additional complications: the integral is 
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not always well-defined. For this reason, it is important first to study such integrals in finite dimension n. 
After this is done, one can take a n — > oo limit. 

However, even in finite dimension n, the averaging operation is not yet fully understood. Most studied are 
Gaussian integrals: that is, when S is quadratic in ^-variables. The simplest Gaussian integral 



(1) = J e *<.'•'•<•'•.'>/";,: 

is easily calculated and expressed through an invariant quantity - determinant - e.g. by diagonalising S: 

1 



J J vAi . . . \ n J 



VdelS 



The integral which factors out is just an 5-independent constant, which can be finite or infinite depending on 
the contour of integration. Therefore, essential S'-dependence of this Gaussian integral is given by (det S)- 1 ' 2 . 
Any other Gaussian integral - with non-homogeneous quadratic S or non-trivial <f> - is equally easy, because 
we have enough freedom to transform quadratic S to diagonal or any other desired form. The possibility of 
diagonalisation greatly simplifies calculations with matrices (tensors with two indices). 

Unfortunately, such methods do not work when S is cubic or higher degree. This can be seen already from 
dimension counting. The number of independent coefficients Syfe, which is n(n + l)(n + 2)/6, exceeds the 
number n 2 — 1 of available SL(n) transformations, so it is not generally possible to diagonalize a cubic action 
- the group SL(n) is too small. That is why the integral 



-s, 



and its higher degree analogues 

Jyi\r j e ^ l! ' n ^ d X, S(xi , . . . , X n ) — Si 1 __i r Xi 1 . . . Xi r (1) 

despite simply- looking, still remain terra incognita. In [1, 2] integral J n \ r (S) was named the integral discrim- 
inant of S, because in the simplest cases this integral is just a power of algebraic discriminant D n \ r (S). For 
example, there is an inspiring formula for 3-forms in two variables 

J 2 ,3 = J e-( ax3+bx2 y +cx y 2+d y 3 )dxdy = (27a 2 d 2 - b 2 c 2 - ISabcd + Aac" + 46 3 d)" 1/6 = (D 2 | 3 ) ~ 1/6 

which shows that J 2 \s and D 2 \3 are, indeed, related. However, when one goes to higher n and r, the relation 
gets more complicated: D n \ r defines only the singularities of J n \ r , while non-singular behaviour is controlled 
by other algebraic invariants. Thus, theory of integral discriminants is closely connected to invariant theory 
[3] and can be viewed as one of the branches of non-linear algebra [2], [4] - [8]. 

It would be very interesting to find a closed formula for generic J n \ r , because it could provide exciting new 
tools in QFT, statistics and other fields where non-Gaussian averaging is used. In this paper, we make a step 
in this direction and find integral discriminants explicitly for 3-forms, 4-forms and 5-forms in two variables, 
and for 3-forms in three variables. Our results indicate an intriguing connection between integral discriminants 
(1) and special functions known as generalized hypergeometric functions [9]. 
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2 Ward identities 



When the action is non-quadratic, diagonalisation and similar linear-algebra tricks fail. To handle non- 
Gaussian integrals, one needs essentially different methods. One of such methods (actually originated in 
the context of quantum field theory) is to find a differential equation, satisfied by the integral as a function of 
its parameters, see [10, 11] for typical applications and references. If such differential equation exists, we call 
it Ward identity (even if it is not directly induced by a change of integration variables). 

Since Ward identities play the central role in present paper, let us give a pair of simple examples to clarify 
this issue. An integral 



m-j 



e - x */2+ax dx 



satisfies a Ward identity 



O -^l F( " ,=0 



because 



(a - F(a) = J (a- x)e~ x2 ' 2+ax dx = J (J- e~ x ' ! / 2+ax ^j dx = 

In contrast with diagonalisation, this method is perfectly generalisable to non-Gaussian integrals: say, 

G(a) = J e- x3 / 3+ax dx 



satisfies a Ward identity 



2 

' G(a) = 



because 



d 2 



G(a) =J(a- x 2 )e~ x3 / 3+ax dx = J (J- e -* 8 /3+a*^ dx = Q 



In this way the problem of non-Gaussian integrals is reduced to another problem - of differential equations. 
This is a much easier problem, especially if differential equations are linear. In the first case we have 



aF(a) = -^F(a) => F(a) = c ■ e"'/ 2 



while in the second case 



d 2 

aG(a) = -Q-pG{a) => G(a) = a ■ Ai (a) + c 2 • Bi (a) 
where Ai and Bi are special functions - Airy functions of the first and second kind. Note that, there is only 
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Arg z — n 



Arg z — 




Figure 1: Two, out of infinitely many, admissible contours of integration for the integral J e x2 / 2+ax dx. 



one linear independent solution in the Gaussian case, while in the non-Gaussian case there are two linear- 
independent solutions. This is because to correctly define an integral, one still needs to specify an integration 
contour. Different integration contours provide different solutions of the Ward identity. 

This relationship between Ward identities and integration contours is quite important, so let us add more 
details. As the simplest option, the contour of integration in the integral 

F{a) = j e- x2 ' 2+ax dx 
c 

can be chosen as C = real axis. However, this choice is by no means unique. In fact, any other contour C 
which asymptotically tends to the lines Arg z — and Arg z = ir/2, is admissible. A contour C is said to be 
admissible, if the integral over C converges. A few admissible contours are shown at Fig. 1. Note that, to 
ensure vanishing of the integral of full derivative (and thus validity of Ward identities) we consider only closed 
contors - the contours at Fig. 1. are closed on the Ricmann sphere, if the infinitely remote point is taken into 
account. Since Ward identity in this case is a differential equation of first order, all the contours give one and 
the same answer F(a) = exp (a 2 /2) up to proportionality. 
More interesting is the non-Gaussian integral 

G(a) = J e- x3 / 3+ax dx 
c 

Note that, the real axis is no longer an admissible contour, since exp ( — x 3 ) grows to infinity when x — > — oo. 
In this case, admissible is any contour C which asymptotically tends to the lines Arg z — 0, Arg z = 2tt/'S 
and Arg z — Att/3. A few admissible contours are shown at Fig. 2. Again, to ensure vanishing of the integral 
of full derivative, we consider only closed contors. Since Ward identity is second order in this case, there are 
two essentially different integration contours, say, C\ and C2. An integral over arbitrary contour - the general 
solution of the Ward identity - is given by linear combination 

J e - x */3+ax dx = Cl J e- x3 / 3+ax dx + c 2 J e- x3/3+ax dx 

Ci Ci 

To summarize the above examples, Ward identity is the main differential equation which governs all the 
contours at once. The choice of particular contour corresponds to the choice of particular solution of the Ward 
identity. For this reason, Ward identities are especially convenient to study of properties, which are invariant 
under change of integration contour. Additional details can be found in [10]. 
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Aig z = — 



Figure 2: Three, out of infinitely many, admissible contours of integration for the integral J e x3 / 3+ax dx. 

3 Integral Discriminants 

Definition. In this paper we study a specific class of non-Gaussian integrals: integral discriminants 

J n \ r (S) = j t " ! ' r/.r | ...,/.r„ 

associated with homogeneous symmetric r-forms S(xi, X2, ■ ■ ■ , x n ). There are two different notations for sym- 
metric forms, which can be useful under different circumstances: tensor notation 

n 

S(Xi, X2 , . . . , ^n) — ^ ^ Si 1 .i 2 ^.^i r Xi 1 Xi 2 . . . X{ r 
i\ l »2,...,» r =l 

and monomial notation 

S{X\ , X2 , • • • j £71 ) ^ ' s a\ ,a 2 ...a„2 ; i 2^2 ■ ■ ■ ^n" 

a 1 +a 2 + ...+a n =r 

To distinguish between these notations, we denote coefficients by capital and small letters, respectively. Note, 
that in tensor notation coefficients have r indices, while in monomial notation they have n indices. 

The choice of contour. Being non-Gaussian integrals, integral discriminants of course depend on the choice 
of integration contour. One has to remember, that not every contour is admissible: for given S, only restricted 
set of contours is allowed. Say, for positive-definite quadratic forms S admissible is any contour, which is 
asymptotically pure real (see the previous section, especially Fig. 1 and Fig. 2 for simple examples). At 
the same time for negative-definite quadratic forms admissible is any contour, which is asymptotically pure 
imaginary. In this paper, we do not describe the contour dependence explicitly. Instead, we concentrate on 
contour-independent properties of integral discriminants. 
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Independence on the form of action. An important feature of integral discriminants, which is due to 
homogeneity of S, is the possibility to substitute the function e~ s under the integral with arbitrary function 
(or, better to say, arbitrary good function) f(S): 



Je ^i.-.^) dx 1 ...dx n ~ J f(jS(x 1 ,...,x n fj dx 1 ...dx n 

Even before specifying the class of good functions, let us consider a simple illustration with f(S) = e~ s2 

f e -(*w) d n x = [ e -(^+-+^l) rl n^^ 1 , f e -(x?+...+x 2 n ) d n x= <-<>nsi 
J J vAiTTTAji J 



(2) 



o -{\ 1 xj + ...+\ n x 2 n ) 



i . . . \ n j 



+•••+< An 



(Fx 



VdeTS 



const 



VdctS 



To specify the class of good functions and prove (2), let us make a change of integration variables 

x\ = p, x 2 = pz 2l x 3 = pz 3 , ...,x n = pz n 

i.e. pass from homogeneous coordinates Xi to non- homogeneous coordinates Zi — Xi/x\. Then 

j e -S(xi,l2 X„) d Xl ...dx n = j P^dp j dz 2 ...dz n e -p r S(hZ2,..,Zn) = 

dz 2 . ■ ■ dz n 



p n - x e- pr d P 



S(l,Z 2 ,...,Z n ) n /r 



For the right hand side of (2) we have 

J f(s(xi,x 2 , . . . ,x n )j dx\...dx n = j p n ^dp j dz 2 ...dz n f(^p r S(l, z 2 , . . . , z n jj 



dz 2 . . . dz n 
S(l,z 2 ,...,z n ) n / r 

Both integrals over p are just S'-independent constants. As one can see, (2) is valid, iff these integrals 

J p n - x e-» r dp and J p^f^dp 

are finite over one and the same contour. This condition specifies the class of good functions f(S). For 
example, all functions f(S) = exp ( — S k ) for k > fall into this class. Relation (2) is therefore proved. As a 
byproduct, we have obtained a non- homogeneous integral representation 



J n \ r (S) = const ■ J — 



dz 2 . . . dz n 



S(l,Z 2 ,...,Z n ) n / r 



(3) 



6 



A useful complement to the definiton of integral discriminants, representation (3) highlights two properties, 
which were less evident in (1): the scaling dimension 



and the vertical symmetry 



J n \ r (\S) — A n l T J n \ r (S) 



In\kr(S k ) ~ Jn\r(S) 



(4) 



(5) 



Ward identities. We now turn to Ward identities, satisfied by integral discriminants with respect to their 
parameters. Let us introduce the correlation functions 

4>{x u x n )^ = j (f>{x l7 Xn ) e -S(xuX2,-,x n ) dxi dXn 



These correlation functions satisfy 



ds a 



■(<t>(xi, . . .,x n f} = -(xl 1 . . .<"0(xi, . . . ,x n )^ 



(6) 



and 



d 



d 



<9s 0l ,...,a n 9s6i,...,6 n 



(j)(xi, . . .,X n )j 



~ai+f>i „,a„+6 r 
j/i . . . 



<j)(xi, . . .,X n )j 



The right hand side depends only on the sum of indices a; + 6j, not on dj and hi separately. For this reason, 
correlation functions satisfy a system of homogeneous second order differential equations: 



d d d d 
dss ds^ dsjj dsg 



4>{xi, . . . ,x n )) = 0, a + b = p + q 



(7) 



In particular, if we set (j){x\ ,x n ) = 1 we find that integral discriminant ( 1 \ = J n \ r satisfies Ward identities 



(8) 



d d d d 



ds s ds^ dsp dsg 



-JJ n \r(S)=0, a + b = p + q 



These equations do not exhaust the set of all Ward identities: there are more. Notice that the differential 
operator in the left hand side of (8) annihilates not only the integral, but also the integrand: 



d d d d 
ds s <9sg dsp ds g 



-S(xi,x a ,...,x„) = Q) a + b = p + q 



and even 



d d d d \ ,/_, ,\ _ ? 

L/Oq L/o ^ fop L/Oq / v 
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which justifies (2) once again. Usually, Ward differential operators do not annihilate the integrand, they just 
transform the integrand into a full derivative (which implies vanishing of the integral over any closed contour) . 
Differential operator in the left hand side of (8) is therefore too special; one can expect other Ward identities 
to exist. To find them, let us consider the most general vanishing correlator - an integral of full derivative 



/ ~Hx~- ' ' ' ' x n)e- S[x " X2 '-' Xn) ) <fci ■ ■ ■ dx n = 



Taking the derivative, we obtain Ward identities for correlation functions: 



/d(j>(x 1 ,...,x n )\ I dS(xi,...,x n ] 



(9) 



For now we are not interested in all correlation functions, only in the integral discriminant. We need to rewrite 
the equations (9) as differential equations on J n \ r . There are several ways to do this, the simplest way is to set 
(j){xi, . . . , x n ) = xj and use the identity (6) to get rid of remaining correlators. Doing so, we obtain equations 



(ySij + flj Aij^J J n \r — 



(10) 



where operators 



■A-ij — ^ ^ s ai...ai~l...aj-\-l...a r , 



d 



ai+(i2 J r...+a n =r 



(ii) 



form the GL{n) algrebra: 



Note also, that other choices of <j){x\, . . . ,x n ) in (9) do not give new Ward identities: everything, that can be 
obtained in this way, will be equivalent to (10). We conclude that the complete set of Ward identities for J„| r 
consists of equations (8) and (10). 

SL(n) invariance. In solving these Ward identities we start from (10), because they are first order. In fact, 
equations (10) simply reflect the Gi(n)-covariance {— SL(n) invariance + correct scaling) of the model. This 
becomes clear, if we separate the generators of GL{n) into two parts: the dilatation (degree) operator 

in + • • • + A nn 



and the other n 2 — 1 operators, which form a representation of SL(n). The first relation implies 

( ^4ll + • • • + A nn J J n \ r = J-n\r 



which is nothing but the scaling dimension (4). The other n 2 — 1 relations imply, that J„i r is SX(n)-invariant 
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function of S. Thus, equations (10) are completely solved by any 5'L(n)-invariant function of S, which is in 
addition homogeneous in S of degree —n/r. A natural question is: is there any simple description of such 
functions? Actually, the answer to this question is positive: any S'L(n)-invariant function can be uniqcly 
represented as a function of the elementary invariants Ik 

Jn\r = F{lk) 

much in the same way as any S , L(n)-invariant function of a matrix Aj can be uniquely represented as a function 
of elementary invariants trA k . Unfortunately, in the case of symmetric tensors S^,. ..,»,. classification of these 
elementary invariants Ik is not that easy, as in the case of matrices. Since a symmetric tensor with r indices 
is much more complicated, than a matrix with two indices, no explicit formula like Ik = trA k is available. 
The study of properties of the elementary invariants Ik, of different kinds of explicit formulas and relations 
between them, is a classical branch of science known as invariant theory [3] . 

The number of SL(n) invariants. Given parameters n and r of the form S, one can easily find the number 
of elementary invariants Ik, of which all other invariants are various functions. Indeed, the linear space of 
forms S of type n\r has dimension 

(n + r-1)! 

dim^i,, = — — - 

r!(n — 1)! 

that is the number of independent coefficients of symmetric tensor S^...^. The group SL(n) acts on this space, 
dividing it into orbits. All forms, connected by SL(n) transformations S^...^ i— > Si 1 ...i r Uj^ . . . ?7j™, U € SL(n) 
belong to one orbit. Ward identities imply, that J n \ r does not depend on coordinates along the orbit ("angular" 
variables). It depends only on transverse coordinates which label orbits ("radial" variables). Simple counting 
of dimensions implies, that dimension of the space of orbits, i.e. the number of radial coordinates, equals to 

i- a j- or/ \ (n + r-1)! 2 

dimi> n | r — aimbL(n) = — n +1 

1 r!(n - 1)! 

Several examples of these numbers are shown at Fig. 3. 



r\n 


2 


3 


4 


5 


6 


7 


2 


1 


1 


1 


1 


1 


1 


3 


1 


2 


5 


11 


21 


36 


4 


2 


7 


20 


46 


91 


162 


5 


3 


13 


41 


102 


217 


414 


6 


4 


20 


69 


186 


427 


876 



Figure 3: The number of functionally independent SL{n) invariants I k of a form of degree r in n variables. 
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The Gaussian case J„| 2 . Note, however, that the case r = 2 (quadratic forms) is exceptional, because the 
dimension dim 12 = n(n + l)/2 is less than the dimension dimSX(n) = n 2 — 1. For this reason, the above 
dimension counting does not work in this (and only in this) case. Actually, as we know, the space of SL(n) 
orbits on quadratic forms is one-dimensional, and the only invariant - the single coordinate on the space of 
orbits - is determinant det S. Thus, determinant is the only variable J n \o can depend on: 

J n , 2 (5) =F(detS) 

The homogeneity condition states, that F(\x) = F(x)/^/\ and has a single solution F{x) = \j\fx. In this 
way we reproduce the well-known Gaussian integral 



l2{s) = vm 



We emphasize, that simplicity of this answer is due to simplicity of the space of orbits, i.e, due to the fact 
there is a single invariant in this case. For higher r, there are many invariants (as many as shown at Fig. 3) 
and J n \ r is a non-trivial function of all of them. The problem is that it becomes impossible to find J n \ r from 
the homogeneity condition alone. The solution of this problem is provided by Ward identities: in addition to 
the homogeneity condition, integral discriminant satisfies relations (8) and this allows to find 

Diagram technique for SL(n) invariants. To do actual computations with SL(n) invariants, we will 
adopt a convenient diagram technique, which is described in [2]. According to [2], it is possible to represent 
tensors with k indices as fc-valent vertices, with two types of indices - covariant and contravariant - represented 
by two types of lines. Contraction of indices is naturally represented as connection of these lines. To simplify 
the diagrams, we use solid lines for both covariant and contravariant indices. In this paper we use the following 
four elementary building blocks for diagrams: covariant tensor 5^...^, represented by a black r-valent vertex 

s 

r indices 

d 

contravariant tensor — , represented by a black r-valent vertex 

r indices 

completely antisymmetric contravariant tensor e* 1 '" 1 ", represented by a white n-valent vertex 



(12) 
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n indices 



and completely antisymmetric covariant tensor e* ± in , represented by a white n-valent vertex 




ri indices 



Tensors e 11 '"*" and in are completely antisymmetric with respect to permutations of their indices, and, 
therefore, remain invariant under SL(n) transformations: 

e il...i nU h yon = ^ U . e h-in = e h-in 

e* i U l > . . . UY = det U ■ e* , = e* , 

where det [/ = 1, because £/ € SL{n). For this reason, any diagram without free (uncontracted) indices, 
made of contravariant e-vertices and covariant S'-vertices, is automatically 5L(n)-invariant function of S 1 . Say, 
determinant ofnxn matrix SV,- can be represented as a diagram at Fig. 4, with n S'-vertices and two e-vertices. 
The other two types of vertices (d/dS and e*) will be used to construct SX(n)-invariant differential operators. 




Figure 4: Determinant of matrix Sij, represented as a diagram of tensor contraction. Black 2-valent vertices 
represent tensor S, white n-valent vertices represent tensor e. 

Diagrams provide a convenient way to obtain SL(n) invariants. In fact, any SX(n)-invariant can be obtained 
in this way, but, unfortunately, not uniquely: absolutely different-looking diagrams can represent one and the 
same invariant. Trying to resolve this ambiguity, one typically faces the complicated problems of classification 
of diagrams and finding relations between diagrams. These problems will not be adressed (and even touched) 
in present paper. Our goal is different: to find an explicit answer for J n \ r in several non-Gaussian cases, using 
diagrams as a convenient tool. We now turn to accomplishing this goal. 
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Figure 5: The degree 4 invariant I4 of a 3-form in 2 variables, represented as a diagram of tensor contraction. 
Black 3-valcnt vertices represent tensor S, white 2-valent vertices represent tensor e. 



3.1 The case of J 2 3 

The simplest non-trivial (i.e. non-Gaussian) example is a 3-form in 2 variables, which can be written as 

S(x, y) = S lllX 3 + 3S 112 x 2 y + 3S 122 xy 2 + S 222 y 2 



By dimension counting, there is only one elementary invariant I4 in this case, given by a diagram at Fig. 5. 
The subscript "4" stands for the degree of this invariant. In this paper we find it convenient to denote the 
elementary invariants of degree k as Ik- It is straightforward to write the algebraic expression for the diagram: 

^4 ^iii2i3 ^313233 ^kik2ka S11I2I3 e ^ C 6 6 6 3 



Evaluating this sum, one gets the following explicit formula for I4 

h = 2S'i 11 S'| 22 — 12S , inS'ii2S'i22'S , 222 + 8S'iiiS , j i 2 2 + 85 , 112 5 , 222 — §S\ 12 S\ 22 

which is nothing but the algebraic discriminant D 2 \ 3 of S: 

D 2]3 = I A (13) 
Since there is only one elementary invariant, the integral discriminant J 2 | 3 must be a function of £> 2 |3 : 

J2\ 3 {S) =F(h) 

Thus, this case is similar to the Gaussian case. The homogeneity condition states, that F(\x) = F(x)/\f\ 
and has a single solution F{x) = 1/ ^fx. In this way we reproduce the formula from the Introduction: 



1/6 

4 



(14) 



In this case, one does not need to use the Ward identities (8): it is enough to use the homogeneity condition. 
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Figure 6: The degree 2 invariant I2 of a 4-form in 2 variables, represented as a diagram of tensor contraction. 
Black 4-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 




Figure 7: The degree 3 invariant 73 of a 4-form in 2 variables, represented as a diagram of tensor contraction. 
Black 4-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 

3.2 The case of J 2 4 

The next-to-simplest example is a 4-form in 2 variables, which can be written as 

S(x, y) = Snnx 4 + 4S 1112 x 3 y + 6S U 22X 2 y 2 + 45' 1222 a;y 3 + £22222/ 

Invariants. By dimension counting, there are two elementary invariants in this case. They have relatively 
low degrees 2 and 3, denoted as I2 and ^3 and given by diagrams at Fig. 6 and Fig. 7, respectively. Looking 
at the diagrams, it is straightforward to write algebraic expressions for I2, 13: 

/., v s' MhMh MhMii 

12 — °JlJ2»3*4 ,J Jlj2j3j4 t t t t 

T„ G .... Q .... C, , Jiji Mh Mki iiki i 3 k 3 jiki 
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Evaluating these sums, one gets the following explicit formulas for l2,h'- 



12 — 2S'iiii5'2222 — 8S'iii25i222 + 6S 112 2 



h — 6S , nii5ii22'S'2222 — 6S'iiiiS' 1 222 — 6S 1112 S 22 22 + 12S'lll2'S'll22'S'l222 — 6iS 112 2 



The algebraic discriminant D 2 \4, just like any other SX(2)-invariant function of S, is a function of I 2 ,h : 

£> 2 |4 = li - 6I£ =8Sf 111^2222 — 965 1111 S'lll25l222S'2222 — 1445 1111 5 1122 S'2222 + 432S' 1111 S'll22'S'i222'5'2222 

— 2165^ 111 5] 1 222 + 432 S mi Sf 112 S 1122 S 2222 ~ 48S , nnS'J 112 S'J 222 5 2 222 — 1440S'iiiiS'iii2S^ 122 <S'i 2 22>S'2222 
+ 6485i ni'S'1122^'2222 + 864S'iiiiS'iii2S'ii22<Si222 + 864S' 1112 S'ii 2 2S'i222'S'2222 + 288>S' 1112 >S 1122 >S' 1222 

— 432S'iiii5'^ 122 S'^ 222 — 216S'] l 112 S'| 22 2 — 512S'j s 112 S'j s 222 — 432Sf n2 S 1 3 1 2 2 £' 2222 (15) 

Integral discriminant. Similarly, the integral discriminant is a function of l2,h'- 

•h\A = F(h,h) 

where the function F is to be determined. The homogeneity condition does not allow to find F unambigously, 
it constrains only the overall scaling factor, but not the non-trivial dependence on the ratio of invariants: 



F(l2,h)=I 2 - 1/A G(± : 



(16) 



Ward identities. To find the function G(z) in this ansatz, we need to use the Ward identities (8). Applied 
to the present case of n = 2, r = 4, the system (8) takes the form 

J214 = 



d 2 


d 2 


ds4Qds22 


ds 3 ids 3 i 


d 2 


d 2 


ds 4a dsi 3 


ds 3 ids 2 2 


d 2 


d 2 


ds4odso4 


ds 3 idsi 3 


d 2 


d 2 


ds 2 2ds 2 2 


ds 3 idsi 3 



J 2 \4 — 



J 2 I4 — 



J214 = 



J2\4 

where s-parameters and 5-parameters are related by 

S40 = Sun, S31 = 4Siii2, S22 = 661122, S13 = 4S1222, S04 = S2222 

Particular equations in this system are, of course, not SL(2) invariant. This is not convenient, since we are 
interested in SX(2)-invariant solutions. Let us transform Ward identities into invariant form, using (16). 
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Invariant form of Ward identities. Substituting the invariant anzatz (16) into the system, we obtain 



d 2 



if 



d 2 d 2 



^2|4 — 



S'll22'S'2222 — "$1222 ((■, . » 2 i a n <9 2 G(z) 1n \^G(z) 



41. 



9/4 



(lUz 2 - 24z) ; ' + (2162 - 12) ■ 



<9s 4 oi9si3 ds 3 ids 2 2 
d 2 8 2 



>Sl222<Sll22 — S'iii2/S'2222 ; . , , ■> 

J?.\i = — • I (l^Az 2 - 24z) 



41. 



9/4 



dz 2 

d 2 G{z) 



dz 



5G(z) 



dz 2 



d 2 G(z) 



(216z-12)^^ + 5G(z) 
dz 

dG(z) 



9s4o9s 4 ds 3 ids 



( d 2 d 2 \ 

\<9s 2 2<9s 2 2 9s 3 i9si3 / 2 



-) J 2 | 4 = ZS ™ S ™~ 3 ^ 122 • ((144, 2 - 24z)^-p -21o: - 1.2, - 5G, :) 



dz 



5iiiiS , 2 222 — S'11126'1222 /,,,,„ 2 ^d 2 G{z) dG{z) 
J 2|4 = TTTgTi I ( 144z - 24*)-^ + (216z - 12) + 5G{z) 



121. 



dz 2 



dz 



where z = I^/li- These equations contain a common SL(2) invariant factor. We conclude, that all the four 
expressions vanish, if and only if G(z) satisfies the differential equation 



d 2 G(z) 



(144z 2 - 24z)^p + (216* - 12) + 5G(z) = 



dG(z) 



dz 2 



(17) 



which is nothing but Gauss hypergeometric equation 

*(! - *) + ic-{a + b+ l)t) abG(t) = 

with a = 1/12, b — 5/12, c = 1/2 and t = 6z. This is the invariant form of n = 2, r = 4 Ward identities. 



Solution. In terms of the Gauss hypergeometric function 

(a + k)T(b + k) T{c) 
T(a) T(b) T(c + k) 



F(\nh\Mt\ f- r(a + fc)r(6 + fc) r(c) t k ab a(a+l)b(b+l) t 2 

^(M.M,*)-!,^ r(6) r( c + fc)M = 1 + T* + c ( c + i) — Y + 



:(c+l) 

the general solution of Gauss hypergeometric equation is given by 

G{t) =ci • 2 Fi([a,6],[c],t)+c 2 -t 1 - c 2 Fi ( [a + 1 - c, 6 + 1 - c] , [2 - c] , i) 



Consequently, the integral discriminant equals 



J 2 , 4 (5) =ci-7 2 - 1/ Vi 





' 1 5 " 




1" 




( 


12' 12 




2_ 


' 4 J 





' 7 11" 




^3~ 


e/ 3 2 \ 


( 


12' 12 




2 





(18) 



where Ci i2 are the two constants, parametrising the general solution of Ward identities. Particular solutions 



4 ) (s)=^ 1/ Vi( 





' 1 


5 " 




"1" 




( 


12' 


12 


1 


2_ 


' 11) 



^) and jW(S) =J 3 / 2 - Y/4 2 Fi 



7/4 





' 7 


11" 




"3" 




( 


12' 


12 




2 


' J2 J 



are associated with different integration contours and can be called the first and the second branches of J 2 | 4 . 
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Singularities. Notice, that the point t = 1 corresponds to 

3 



which is just the discriminant's vanishing locus if — 6if = Z? 2 |4 = 0. This is interesting, because the point 
t = 1 is a singular point of the hypergeometric function 2 -Fi- However, there are two other suspicious points: 
hypergcometric function 2-F1 can have singularities at t = 0, 1 and 00. Let us study asymptotics of integral 
discriminants at these points, using the formulas 

2 F 1 ([a,b],[c],t) = l + 0(t), whent^O (19) 
2 F, ([a, b) , [a + b] , 1 - t) = - r ( ( °)r(b) 1 °S t + °( t0 ) ' wh ent^0 (20) 
2*1 ([a, 6], [c], t) - f ~" rjb)r(t-a) ' 2i?1 ([ °' - c + !]' [° - 6 + !]. V*) + 

^ r(a)r(c-fr) ' 2Fl 6 - c + 1], [6 - a + 1], 1/t) (21) 

which can be found in any reference book of hypergeometric functions (see, e.g. [12]). Using (19), we find 

J 2|i ( S ) ~ *2 V4 + ° (*) . when * -» 
J 2|i ( 5 ) ~ V4 + O(t) , when t -» 

Using (20), we find 

4' W ~ r(i7i r 2) 1 rwi2) J ' 1/< 108 (1 - f > + o((0> - when ' - 1 

■*S (s) ~ mmwrn/^'' los (1 - f > + o( ' 0) • when « - 1 

when i — > 00 
when t — > 00 

We can see from these results, that singularities of the integral discriminant J 2 | 4 at finite values of h,h are 
completely controlled by the algebraic discriminant £> 2 | 4 . The singularity at discriminant locus is logarithmic. 
Other singularities, not related to the algebraic discriminant, are situated at infinite values of invariants: 
I2 — > 00 (t — ► 0) and J3 — > 00 (t — ► 00). 



Using (19) and (21), we find 



4?AS) - M (-6)-^/ 3 - 1/6 + Q(r 5 / 12 ), 

^2(s) - TS 2) (- 6 )" 7/12j 3 1/6+Q ( rll/12 )' 
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Hypergeometric integral. Function 2 F\ has an integral representation 



2 F 1 ([a,b},[c\,t) 



F(c) 



T(b)T(c 

which is a direct consequence of a simpler identity 



s^il- s) c - b - 1 



ds 



(i - sty 



(22) 



(23) 



For the integral discriminant, we obtain 



'2|4 



1 s -7/12 (1 _ s) -ll/12 i s -l/12 (1 _ s) -5/12 

as — — — h C2/3 • / «s 



(Jf - 6s/I)Vi2 



(24) 



Hypergeometric functions with different values of parameters a, b and c are related by various "modular" 
transformations of the variable t and the integration variable s, which leave the boundary region (two points 
and 1) intact. Such transformations were first found and studied by Euler, therefore they are known as Euler 
hypergeometric transformations. Say, transformation 1 1— > t, sh^I - s gives rise to a relation 



2 F 1 ([a, b] , [c] , t) = (1 - t)- a ■ 2 F 1 ([a, c-b],[c\, j^j 



(25) 





' 1 


1 " 




1" 




( 


12' 


12 




_2_ 


' D 2l J 



"7 7 " 




3" 




12' 12 




_2_ 


' ^2! 4 ; 



After this transformation, J 2 | 4 takes form 

J$(s) = (D 2l4 y 1/12 2Fl 
j$(s)=i 3 (D 2li y 7/12 2 F 1 

In this way the hypergeometric integral (24) allows to recast J 2 |4 m various forms and establish relations 
between them. Therefore, (24) is a useful and important representation. We emphasise, that the relation of 
integral J e~ s ^ x,v ^dxdy to hypergeometric integrals (24) is not a priori expected: these integrals look very 
different, even from the point of view of variables they depend on. 

Vertical symmetry. To finish this section, let us check the vertical symmetry between J 2 | 2 and J 2 | 4 : 

J 2 | 4 {{ax 2 + bxy + cy 2 ) 2 ^j ~ J 2 | 2 (ax 2 + bxy + cy 2 ^j 

Setting S(x, y) = (ax 2 + bxy + cy 2 ) 2 we obtain the coefficients 



Sim = a 2 , S m 2 = -ab, Sn 22 = -ac + -b 2 , Si 222 = -be, S 2222 = c 2 
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and the invariants 



h = \{b 2 -±acf, / 3 = l(62_ 4ac) 3 

00 



Substituting them into (18), we find 



where the constants of proportionality 
A (i) = 6 -i/4 . 2Fi 

A (2) = g+1/4 . ^ 



' 1 5 




1 — 1 


12' 12 


; 


2 


" 7 11" 




3" 


12' 12 


5 


2 



Vfo 2 - 4ac 



1 



r(i/2) 



r(i/i2)r(5/i2) 



.6+1/4. 



r(3/2) 



log <~ oo 



■ log <~ oo 



r(7/i2)r(ii/i2) 

are independent of a,b,c- as prescribed by the vertical symmetry - but infinite. Only their linear combination 

finite number = — — 



(Vim r(3/2) \ (1) _ / 1/4 r(i/2) \ (2) 
^ r(7/i2)r(ii/i2); ^ r(i/i2)r(5/i2); 



is finite logarithmic divergencies cancel out. Actually, it is expectable that only one linear combination of 
solutions stays finite. This is because not every contour, admissible for J 2 |4, is admissible for J 2 | 2 . Admissible 
contours for J 2 \4 approach infinity from 4 different directions, where the 4-form (ax 2 + bxy + cy 2 ) 2 takes real 
positive values. Only two of these directions are appropriate for J 2 | 2 - those, where the 2-form ax 2 + bxy + cy 2 
takes real positive values. One linear combination of branches corresponds to admissible contours, while the 
orthogonal linear combination corresponds to other contours. That is why both branches are singular at 
t — ► 1, but certain linear combination of branches is regular at t = 1. Note, that nothing similar happens at 
non-discriminantal singularities: at t — > and t — > oo no linear combination of branches stays regular. 

3.3 The case of J 2 \5 

Our third example is a 5-form in 2 variables, which can be written as 

S(x, y) = Smiix 5 + 5S*i m2 a; 4 y + 10S , m 22 x 3 y 2 + 10S*ii 222 a; 2 y 3 + 55i 2222 xy 4 + 5 22222 y 5 

Invariants. By dimension counting, there are three elementary invariants in this case. They have degrees 
4, 8 and 12, denoted as I4, Is and ii 2 and given by diagrams at Fig. 8, Fig. 9 and Fig. 10. Looking at the 
diagrams, it is straightforward to write an expression for I4 



T. — Q Q <?,,,,, <?,,,,, J1hJ2j2J3j3J4.k4.J5H jili j 5 h k-Lh k 2 l2 f k 3 l3 

± 4 u lll2l3i4i5 u JlJ2J3J4Jr>klk2k3k4k' ) u Llt2l3Uh C CCCCCCCCC 



and equally straightforward to write expressions for J 8 ,/ 12 . Evaluating the contraction, one gets a formula 
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Figure 8: The degree 4 invariant I4 of a 5-form in 2 variables, represented as a diagram of tensor contraction. 
Black 5-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



I4 = 2S , 11111 S , 22 222 — 20SiiiiiS , iiii2<Sl2222<S'22222 + 8£lllll'S'lll22'S'll222<S22222 + 

325 , iiiiiS , 11122<S'i2222 — 245 , iiiii5'i 12 22<S'l2222 + 325 , j 1112 5 , 11222<S'22222 + 185 , iiii 2 5'l2222 — 

24S , nii25'i 11 22'S'22222 — 152S'llll2/S'lll22<S'll222<S'l2222 + 96SlH12<S'ii222 + 965 , 111 2 2 S'l2222 — 64S l i 11 2 2 S'ii222 

and similar formulas for Is,Ii2 - they are quite lengthy and we do not present them here. The algebraic 
discriminant D 2 \5, just like any other S'L(2)-invariant function of S, is a function of I4,l$,li2'- 

D 2 |5 = II - 64/ 8 (26) 
Integral discriminant. Similarly, the integral discriminant is a function of Ii,I%,I\2'- 

J2\5 = F(h,I 8 ,I 12 ) 

where the function F is to be determined. The homogeneity condition does not allow to find F unambigously, 
it constrains only the overall scaling factor, but not the non-trivial dependence on the ratios of invariants: 

^ 4 ,/8,/ 1 2)=/ 4 - 1/10 G(|,^) (27) 
To find the function G(u, v) in this ansatz, we need to use the Ward identities (8). 
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Figure 9: The degree 8 invariant Is of a 5-form in 2 variables, represented as a diagram of tensor contraction. 
Black 5-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



Ward identities. Applied to the present case of n — 2, r — 5, the system (8) takes the form 



d 2 


d 2 


ds 50 ds 32 


ds4ids4i 


d 2 


d 2 


ds 50 ds 2 3 


ds4ids 32 


d 2 


d 2 


ds 50 ds u 


ds 41 ds 2 3 



J 2 \5 = 
J 2 \5 = 
J 2 \5 = 



d 2 


d 2 


ds3 2 d.S32 


ds4ids 2 3 


d 2 


d 2 


ds5odso5 


ds 32 ds23 


d 2 


d 2 


ds4idsi4 


ds 32 ds23 



where s-parameters and 5-parameters are related by 

s 50 = Sum, S41 = 55mi2, S 32 = IOS11122, «23 = 10Sll222, «14 = 5S'i2222, S05 = 5*22222 



Particular equations in this system are, of course, not SL(2) invariant. 



Invariant form of the Ward identities. It is possible to deal with non-invariant equations (as we did in the 

previous section) but it is much more convenient to rewrite the above differential equations in SX(2)-invariant 
form, by summing them with appropriate 5-dependent coefficients: 



Oo = (2S'iiiiiS'i2222 — 8S11112S11222 + 6S11122) 2— — — 8— — — h 6- 



ds^o dsi4 ds4i ds23 ds 32 <9s 32 
(33 3 3 3 3 

2(S'iiiiiS , 22222 — 3S , ini2/S'i2222 + 2SiH22<S'll222) TJ TJ 3— — h 2- 



ds 05 3s 5Q 3s4i 3si4 3s 32 ds 23 
,3 3 3 3 3 3 

(2S'iiii 2 S' 2 2 2 22 — 8S'iii22S'l2222 + 66^222) ( 2 — — 8— — h 6- 



i9s 4 i 3s 05 3s 3 2 3si4 3s 2 3 ds 23 
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Figure 10: The degree 12 invariant J12 of a 5-form in 2 variables, represented as a diagram of tensor contraction. 
Black 5-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



and 



4 =Pll(S) 




d 




a 




ds 14 


ds 4 i 


<9s 2 3 


2Pi 2 (5) 


( 9 


d 


-3/ 


a 


\ds Q5 


ds 50 


dsn 


<9si4 


P22(S) 


f 2 9 


a 




a 


\ ds 4 i 


ds Q5 




5si4 



,5 a , 

+ 6^ "7: + 



„ a a . 

+ 2^ "7: + 



6-t — 



where P is a quadratic form with coefficients 

Pab(S) — S a i 2 i 3 i 4 i 5 Sj 1 j 2 j 3 j 4 j 5 Sk 1 k2k 3 k4k^Sl 1 l 2 l 3 l i l 5 S mim2rn3m4m5 S SlS2S3S4 l, 



e *2j2 £ i3i3 £ ufc4 e i5k 5 e kih e k 2 h e j4m4 e j5m 5 e miS! e m,2S2 e hs3 e hs4 e jik 3 e l 5 m 3 



Operators Oo and O4 are SX(2)-invariant, simply because they are given by diagrams at Fig. 11 and Fig. 12. 
The subscripts "0" and "4" stand for the degrees of these operators. In analogy with invariants Ik, we denote 
the invariant differential operators of degree k as Ok- 
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Figure 11: Invariant differential operator Oo for the case 2|5, represented as a diagram of tensor contraction. 
Black 5-valent circles represent tensor S, white 2-valent circles represent tensor e, black 5-valent rhombuses 
represent tensor d/dS, white 2-valent rhombuses represent tensor e*. 




Figure 12: Invariant differential operator O4 for the case 2|5, represented as a diagram of tensor contraction. 
Black 5-valent circles represent tensor S, white 2-valent circles represent tensor e, black 5-valent rhombuses 
represent tensor d/dS 7 white 2-valent rhombuses represent tensor e*. 
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Of course, operators Oo and O4 are not unique - there are many other invariant differential operators, which 
annihilate J 2 |5- All such operators contain an e-antisymmetrized combination of two derivatives d/dS, which 
is exactly the structure of Ward identities. We do not study here the interesting problems of classification of 
these operators and finding relations between them, because Oo and O4 are quite enough to find J215. Let us 
derive a useful formula for the action of second order operators 



a,f} 



_d d_ 



By application of the chain rule, we obtain 

k k a,f3 



dh 'ri ' p dS a dS fj f-' dl k dl, 



\ - d 2 F x - dl k dl m 



k,m 



ct,0 



dS a dS 



Since 



a,/3 



d 2 h 

8S a dSp 



Oh 



and 



^ 2C Q ,,^— 7T-TT FT = O(hlm) - IkOI m - I m OIk 
a.fj 



dl k dl, 



we finally obtain an important formula: 



dF 



d 2 F 



2 dl k dl„ 

k.m 



6{l k I m ) - IkOIm - ImOh 



(28) 



As one can see, to describe the action of O-operators on arbitrary invariant functions it suffices to calculate 
the action on all invariants h and all products hl m - In the present case we have 



/ \ 

h 

h 

1X2 
\ ) 



( 



264 



\ 



25 h 



12 
25 



294 
~25~ 

162 



4-<8 



'12 



74 7 8 7 4 7i 2 



I&Il2 7 12 



25 



384 . 
~25~ 



2 , 1166 192 

— 7? H 7 4 7 8 H A 

25 4 25 5 



12 , 144 , 1794 
-ll h + —lt + —UI 12 



25 



1166 192 



8 , 1188 , 12 

-/I/8 + — /i + -/4/ 12 



9 2 12 2 1944 

^^4/ 8 +-/4A2 + — /8/12 



12 , 144 , 1794 



9 , 12 , 1944 

-7 4 7| + -,fl 12 + _, 8 / 12 



54 , 84 684 , 

-_/| + _/ 4 7 8 7 12 + _/ ? 
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and similarly for the second operator 



4 



/ \ 

h 

h 

hi 
\ ) 



264 

2 r r 588 t 

363, 2 153 
~W h ~ "25 ' hh2 



64 



J 4 


/4/ 8 


/4/12 




A 2 




/4/12 




I 2 
1 12 



V 



928 768 
~125~ 4 8 ~ 1>5~ 12 


2 , 
"25^ 


584 2 ,524 
" + ~25 hh2 


363 2 

^T /4/s ■ 


153 2 
" "^5 4 


696 
/l2 " IT* 


2 2 584 2 524 
"2^ /478 "l5" /s + l5" /4/l2 


— / 4 /| - 
25 8 


22 2 2376 
2^ /4/l2 + ~ /8/l2 


603 /3 
50 8 


291 

Uh 

25 


1188 
/12 + — 


363 2 153 f2 , 696 


603 /3 
50 8 


291 1188 2 
— /4/8/12 + — 1 12 


129 
~5~ 


^12 - 


606 2 
^5" 74/12 



Applying (28), we obtain 
O F(. 



'2 2 294 \ <9F /12 162 \ 9F 



+ ( 8 ^ 4 2 



192 
~25~ 



d 2 F 

~dif 



|/l/ 8 + m| + ^i 2 )g 



27 
— j 

25 



,608 T r 192 



d 2 F /144 9 144 \ a 2 .? 1 



3 

25" 



+ ( -T^hlj + 7/^12 + ^hh2 



168 
T 



d 2 f 
dhdh 2 



and similarly for the second operator 



0^,4,^) = -f/ 8 g + (. 



2 588 T \ <9F /363 r2 153 r r 

— hh + —hi ^ + — II - —hi 



OF 

25 ^ 12 J^ 



384 \ d 2 F ( 1 
'12 ) 7T72 



+ ( -87 , /, - —l yi ] — + j ^7 4 7| - ilj|j 12 + 24/ 8 7 12 



25 



97? 



10 



25 



^T 1 /141 



97 2 + V 25 



■^8^12 — 6/ 4 7 2 2 



d 2 f 

U1 12 



+ 



64 _ 64 
" 5 h 25 



. , , d 2 F 432 <9 2 F 
^4Ji2 nr ar TrrHhi- 



dhdl s 25 ^d7 4 <97 



12 



24 A 3 136 AAA + 94A 2 



d 2 F 

dl 8 dl 12 



For the integral discriminant, both expressions vanish. Substituting the ansatz (27) and making the necessary 
algebraic transformations, we obtain two differential equations on the function G(u,v): the first 



d 2 G d 2 G 
50(-l + 64m)(m + 6m 2 + 15m)— r- + (75u 2 + 72000m 2 + 57600mm 2 + 600mm + 7200m 3 - 250m)—— + 

ou l ouov 

d 2 G dG 
(675m 3 - 13500w 2 + 10800mm 2 - 750mm + 43200mm 2 )—— + (50400m + 30720m 2 - 50 + 5770m) — + 

ov z ou 

dC 

(-300m + 60480mm + 11160u 2 - 7650m)— + (528m + 110)G = 
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and the second 

d 2 G d 2 C 
25(-l + 64w)(5w 2 + A8uv - 22v) — T + (230400m; 2 - 39600i; 2 - 6000u 3 + 28800uu 2 + 6800wm)— - + 

d 2 G dG 
(172800w 3 + 750Cto 2 + 25200m; 2 - 7050vu 2 )—^- + (-36120w + 4000u 2 + 122880im + 100u)— + 

dG 

(241920w 2 + 19440wu - 9075w 2 + 7650i;) — + (-220u + 2112u)G(u, v) = 

where u = I$/I 2 and v = These two linear differential equations in two variables constitute the 

invariant form of n = 2,r = 5 Ward identities. Integral discriminant is found as the solution of this system. 



Solution. Having linear differential equations, it is easy to solve them in series: if one puts 



G(u,v) = ^ CjjU l yi 



then one finds, after some algebraic transformations, 
Ca = const • — — • 16 



r(- + i + 2ijr(- + i + 2i 

i.e. there is unique series solution. In this way one obtains the first branch of the integral discriminant J 2 \5'- 




(29) 



This answer is interesting: the function appears to be of generalised hypergeometric type [9] , since its coeffi- 
cients Cij are ratios of T-functions, depending on linear combinations of indices i and j. However, these series 
have a rather small convergence radius. Already the subsequence with i = j behaves as 




and diverges when the combination in the last brackets exceeds unity. For analytical continuation, one better 
substitutes this series by its integral representation. 



Hypergeometric integral. Applying the formula (23) twice, we get 



10 



r TFT + • + 3 r + 2i + 3 i r +3 



10 



10 



r - + i + 2j T - + i + 2j 



G 



1 1 

// 



dtds 



t-"" W {8-t8)-»' W 



t(l-t)s c 



25 



The sum over i and j is calculated, using 



hj r i - - 1 - j 



A 1 B J _ 1 

X7T7T ~ vi + a + b 



and we obtain the following integral representation: 



1 1 

G(u, v) = J J dtds 





t -7/10 s -9/10 (1 _ t j-9/10 



V3 - 3s + A8uts 2 + 128vsH - 128i;s 3 t 2 



(30) 



Accordingly, for J 2 i 5 we have 



l l 



t -7/10 s -9/10( 1 _ t )-9/l 



^3/1 - 3/|s + 48hl$ts 2 + 128I 12 sH - 128/i 2 s 3 t 2 



(31) 



Other branches of J 2 |5 can be obtained by various transformations of variables (u, v) and integration variables 
(t, s), which leave the boundary region intact ("modular" transformations). We do not list them in this paper. 

Singularities. One of the benefits of integral representation is the possibility to analyse the singularities. 
Singularities of the integral (30) come from zeroes of the polynomial 



P(s, t) = 3 - 3s + A8uts 2 + 128vsh - 128vsH 2 



which stands in the denominator. Ordinary zeroes (say, in t) are inessential, already because the singularity 
dt/y/t — to is integrable, and - more important - because the integration contour can be deformed and taken 
away from the singularity. This can not be done if the two roots coincide, i.e. if discriminant of P(t) is 
vanishing. Then what matters is integration over s. If 

Disct(P(s,t)) - (t^s) -t 2 {s)f 



has a simple zero at s = s , then 

P(a, t)~(t- t Q ) 2 + (h - t 2 ) 2 ~ (t - t ) 2 + (a - s ) 



and the integration contour for s can be taken away from the singularity. More serious is the case when 
discriminant has double zero, i.e. when repeated discriminant vanishes, 

Disc s ( Disct(P(s,t))) = 
Then contour can not be deformed, neither in t nor in s. However, in the vicinity of such point 
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Figure 13: The degree 4 invariant I4 of a 3-form in 3 variables, represented as a diagram of tensor contraction. 
Black 3-valent vertices represent tensor S, white 3-valent vertices represent tensor e. 



P(s,t) ~ (t - t ) 2 + (s - s ) 2 



and this singularity, though unavoidable, is integrable because of the square root and because the integral is 
two-dimensional. This is most simply expressed in polar coordinates: 

dsdt 

drd<p 



Moreover, even if repeated discriminant has zeroes ofhigher order, P(s,t) <~ (t — 1 ) 2 + (s — s n ) k , the singularity 
in the integral does not arise. 

Thus the only remaining source of singularity is when it occurs at the ends of integration segments (at 
the boundary of integration domain - because the contour can not be deformed to avoid these boundaries), 
i.e. when s or t equals or 1. Given the shape of P(s, t), of these the only significant one is at t = 1, when 
P(s, 1) = 3 — 3.s + 48us 2 . Remaining integral over s becomes singular when this expression has a double zero: 

Disc s (3 - 3s + 48hs 2 ) = 1 - 64u = 
i.e. when u = 1/64. This point corresponds to 

64J 8 _ 

r 2 

which is just the discriminant's vanishing locus J| — 64/g = D 215 = 0. We encounter once again the same 
relation between integral and algebraic discriminants: the latter controls essential singularities of the former. 
Note, that the double zero occurs at the point s = 2, which lies beyond the integration domain. For this 
reason, the singularity is rather soft. 



27 



Figure 14: The degree 6 invariant Iq of a 3-form in 3 variables, represented as a diagram of tensor contraction. 
Black 3-valent vertices represent tensor S, white 3-valent vertices represent tensor e. 

3.4 The case of J 3 3 

Our final example is a 3-form in 3 variables, which can be written as 

S(x, y) = Siiix 3 +3Sii2x 2 y+3Sii 3 x 2 z+S222y 3 +SSi22xy 2 +3S223y 2 z+S 3 3 3 z 3 +3Si3 3 xz 2 +3S233yz 2 +6Si23xyz 

Invariants. By dimension counting, there are two elementary invariants in this case. They have degrees 4 and 
6, denoted as I4,Iq and given by diagrams at Fig. 13, Fig. 14. Looking at the diagrams, it is straightforward 
to write the algebraic expressions for 7 4 and I e : 

d 4 — ^Z 1 l2l3 U Jlj233 U Klk2k3 tJ lll2h C tec 

r„ — C?. ..<?...<?,, , <?,,,<? <? Jikili c i232S2 c 3ik2rn 1 l 2 rn 2 k3 m 3 s 3 j 3 l 3 i 3 s 1 

J 6 iJ l\l2l3 U 3l32j3 u KlK2 K3 ^ J l-\_l2l3 iJ m\m2m3 tJ S\S2S3 Ki t c t t t 

Evaluating these sums, one gets the following explicit formulas for 

I4 = 6Sf 2 3 — 12S'i22<S'i23>S'l33 + 65'i22>5'l33 + 6Su 3 Si23Si 33 S222 — ^2Sn 3 S 2 23 S223 — 6Sii3>S'i22<Sl33<S'223 + 
18S'ii35i225'l23'S'233 — 6 ( Sii3 ( Sf 2 2'S'333 + ^S 2 13 S 2 23 — 6Sf 13 S'222<S'233 ~ GSu2S 2 33 S222 + 18S'll2<S'l23<S'l33<S223 — 
12S'ii25 , j23'S'233 — 6S'ii2S'l225'l33'S'233 + 65'ii2S'l225'l23'S'333 — 6S'ii2iS'll3iS'223>S , 233 + 6S'ii2S'll3iS'222<S , 333 + 
6S'ii2'S'233 — 6S'i 1 2'S'223'S'333 — 6S'iii/S'i33/S'223 + 6S'iii/S'i33/S'222'S'233 + 6S'iiiS'i23'S'223'S , 233 — 6S'iii/S'i23<S'222<S'333 — 

65iiiS'i22<S , 233 + 65'iiiS'i22'S'2235 , 333 

la = 48^23 — 144S'i22<S'i23'S'i33 + 144Sf ^S^gS^g — 485*^22 S 133 + 72Sn 3 S 3 23 Si 33 S222 — 144Sii3S ,4 2 3>S , 223 — 

72S , ii35i22<5i23S'i33S'222 + 72Sli3<Sl22<S'i23<S'l33<S'223 + 216S'ii3»S , i22<S , i23<S'233 + 72Sii35f 2 2'S'l33'S'223 — 

216S , ii3S'i22'S , i235'i335'233 — 72S\ 13 S 2 22 Sj 23 S333 + 72Sn 3 S 3 22 Si 33 S 333 + 18S 2 13 S 2 33 S 2 22 - 
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72Sf 13 <S'i23£'i33£'222£'223 + ^^^113^123^223 — 72Sf 13 Sf 2 3S222S 2 33 + 725i 13 5i225i335f 2 3 - 
36S , ii3»S'l22'S'l33'S'222'S , 233 — 216S , ^ 13 S'i22S'l23'S'223'S'233 + 1445i 13 5i225i23 5222 5333 + 162S'j 13 S'j22'S , 233 — 
1445 f ^ 13 S , f 2 2'S'223'S'333 ~ 48S'j i 13 S'| 23 + 72S'j i 13 S , 222<S'223'S , 233 — 245 3 13 5f 2 2 $333 — 72S'ii2iSf 23 ,Sf 33 S'222 + 
2165'll25 , f 23 S'l33S'223 ^ 144S'll25 , ] 1 23 S'233 + 72£ll2<S'l22<S'i33<S , 222 — 2165ll25l225l23 5i 33 5223 + 
72S' 11 2>Si22<S'i23>S'l33'S'233 + 72S , 112 5'i225'i23'5'333 + 72 5 f 112 S'f 22 S'^ 33 S'233 — 72S\ 12 Sf 2 2 $123 ^133 5333 ~ 
365ii25ii3S'j 33 5'222<5 , 223 — 216S'ii2£'ll3<S'l23<S'l33<S'223 + 360S'ii25ii3S'i23<S'i335'222'S'233 + 72S'ii25ii35^ 23 5'2235 , 233 — 

2165 112 5 

113'5l23'5'222'S'333 + 36S , ii 2 5'ii3S'i22'S'l33'S'223'S , 233 — 108S 112 S 113'S'l22'S'l33'S'222'5'333 — 
216S'ii2S'ii35i22'S'l23'5'233 + 3606'ii2S'll3/S'i225'l23'S'223'5333 — 36S'ii25ii3S'^ 22 S' 2 33S , 333 + 72S'ii2S'i 13 S'|2 3 5'233 — 
144S'ii2S'i 13 S'222'S'233 + 72S'ii25'i 13 S'222'S'223'S'333 + 1625 112 5 133 S 22 3 ~ 1445ii 2 5i 33 5222 5233 — 
2165'J 12 5'i23'S'l33'S'223'S'233 + 144S'^ 12 5'i235'l33'S'222'5'333 + 144S , i 12 S'i 23 5| 33 — 72Sf i2<5'i23'S'223'S'333 + 
72S^ 12 S'i22<S'l33>S233 — 36S'j 12 S'i22/S'i335'223'S'333 — 72Sf 12 <S'l22'S'l235'233'S'333 + 185^ 12 Sf 22 5f 33 + 
725i 12 5ii35223'S'233 — 144S'i 12 5'ii35'| 23 5333 + 72S'f 12 5'ii3S'2225'233'S'333 — 48S , f 12 S' 233 + 725 , f 12 S , 223'S'233'S'333 ~ 
245'ii2'5'222'5'333 ~ 245*11151335122 + 1445iii5i235i 33 52225223 — 725iii5i 23 5i335|23 — 
2165in5^235i3352225233 + 725m5 3 23 5223<S'233 + 1205in5^2352225333 — 1445iii5i225^3 3 5| 2 3 + 
725iii5i225i 33 52225233 + 3605in5i225i235i3352235233 — 725m5i225i235i3352225333 — 725lll5l225 12 35 23 3 — 
216511151225^2352235333 — 1445m 5^2 >S'l33>S233 + 725m 5^2 >S'l33>S'223>S'333 + 1445in5^ 2 2'S'l23'S'233'S'333 — 
245iii5i225| 33 + 725in5ii35i335|23 — 1085m5n35i335222'S l 223'S l 233 + 365iii5ii 3 5i 33 5f22'S'333 — 
725ni5ii 3 5i235| 2 35233 + 1445iii5ii35i23 5222 5f 33 — 725i Il5n 3 5i23522252235333 — 365m 5ll35i22 5223 5| 33 + 
725m5n35i225| 2 35333 — 365iii5ii35i22 5 2 22 52335'333 — 365m5ii25i335|2 3 5233 + 725m5ii25i3352225| 3 3 — 

365m5ii25i33522252235333 — 725m5ll25l2352235 2 3 3 + 1445m5ii25i235| 2 3 5333 — 
725m5ii25i23522252335333 + 725m5ii25i225f 33 — 1085m5ii25i22522352335333 + 365m5ii25i2252225| 3 3 + 
185m5 22 35 2 33 — 245m5 2 2 3 5333 — 245 m 52225 23 3 + 365 m 5222522352335333 — 65 m 'S'|22 , 5'333 

The algebraic discriminant -D 3 | 3 , just like any other 5L(3)-invariant function of 5, is a function of J 4 and I 6 : 

D 3{3 = 32/| + 3/ 2 (32) 

When expanded, discriminant -D 3 | 3 contains 2040 monomials. See the appendix of the book version of [2], 
where it is written explicitly for curiosity. Formula (32) is a remarkably concise expression of this disriminant 
through a pair of invariants, given by beautiful diagrams Fig. 13 and Fig. 14. It is interesting to extend this 
type of formulas - eqs. (13), (15), (26) and (32) - to higher n and r. 

Integral discriminant. Similarly, the integral discriminant is a function of I^.Iq: 

J m = F(h,I 6 ) 

where the function F is to be determined. The homogeneity condition does not allow to find F unambigously, 
it constrains only the overall scaling factor, but not the non-trivial dependence on the ratio of invariants: 

F(h,I 6 )=I- 1/4 G^ (33) 
To find the function G(z) in this ansatz, we need to use the Ward identities (8). 
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Ward identities. Applied to the present case of n — 3, r = 3, the system (8) takes the form 
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where s-parameters and 5-parameters are related by 

•S300 = Sin, S210 = 35*112, S201 = 35*113, s 030 = S'222, Sl20 = 3Si22, S02I = 35*223 
s 003 = 5333, S 10 2 = 35i33, SQ12 = 35233, s lll = 65i23 



Particular equations in this system are, of course, not SL(3) invariant. 

Invariant form of Ward identities. To rewrite the above differential equations in 5L(3)-invariant form, 
we use again the method of invariant differential operators - sum the Ward operators with appropriate 5- 
dependent coefficients to form an invariant operator O4, given by the diagram at Fig. 15. The antisymmctrized 
combination of two derivatives (two black rhombuses) which is a part of this diagram, ensures that O4 is a 




Figure 15: Invariant differential operator O4 for the case 3|3, represented as a diagram of tensor contraction. 
Black 3-valent vertices represent tensor 5, white 3-valent vertices represent tensor e, black 3-valent rhombuses 
represent tensor d/dS, white 3-valent rhombuses represent tensor e*. 
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linear combination of Ward operators and annihilates J 3 | 3 . Operator O4 contains 4 horizontal lines in the 
inner circle and belongs to an infinite family of operators with 2p horizontal lines in the inner circle: 



2p , p= 1,2,3,... 

of which the simplest are Oq and O2, given by diagrams at Fig. 16. However, operators Oq and O2 are too 
simple: in fact, they do not constrain J 3 | 3 at all. Thus O4 is the main operator we use in this section. 





Figure 16: Invariant differential operators Oo and O2 for the case 3 1 3, represented as diagrams of tensor 
contraction. Black 3-valent vertices represent tensor S, white 3-valent vertices represent tensor e, black 3-valent 
rhombuses represent tensor d/dS, white 3-valent rhombuses represent tensor e*. These operators cannot be 
used to find J 3 | 3 , because they do not give any non-trivial equation on it. 



The action of O4 on invariants and their products is given by 
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Applying (28), we obtain 
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For the integral discriminant, this expression vanishes. Substituting the ansatz (33) and making the necessary 
algebraic transformations, we obtain a differential equation on the function G(z): 



(144z 2 + 1536z) ^ ( 2 Z) + (216* + 768) + 5G(z) = 



dz 



(34) 



which is nothing but Gauss hypergeometric equation 



*(! - + (c-(a + b+ l)t) ^ - abG(t) = 



with a = 1/12, b = 5/12, c = 1/2 and t = — 3/32z. This is the invariant form of n — 3, r = 3 Ward identities. 
Solution. In terms of the Gauss hypergeometric function, the integral discriminant equals 



J 3 \ 3 {s) = Cl -r 4 1, \f 1 
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(35) 



where c lj2 are the two constants, parametrising the general solution of Ward identities. Particular solutions 
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and J<j$(S) 
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12' 12 
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' 32/ 4 3 J 



are associated with different integration contours and can be called the first and the second branches of J 3 | 3 . 
One can see, that the results of this section are parallel to those of sec. 3. 2 - even the rational parameters in 
the hypergeometric function are the same. 

Singularities. Notice, that the point t = 1 corresponds to 

'Ml . 



32/ 4 3 



which is just the discriminant's vanishing locus 32/1 + 3/g = D 3 | 3 = 0. To investigate the two other suspicious 
points t = and t = oo, let us study asymptotics at these points. Using (19), we find 

J^l (S) ~ J7/ 1/4 + O (t) , when t -» 
(S) ~ I 6 i; 7/i + 0(t), when t -» 



Using (20), we find 



J » s ) - iwiiW" ,4log<1 + mi* +0<t ° ) ' wte ° ' "* 1 
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Using (19) and (21), we find 

W~ J ^(£rV'- + o(^>. w he „ 

^^^^(^'"^^n when. 



t — > oo 



oo 



We can see from these results, that singularities of the integral discriminant J 3 | 3 at finite values of Zi, Jg are 
completely controlled by the algebraic discriminant -D 3 | 3 . The singularity at discriminant locus is logarithmic. 
Other singularities, not related to the algebraic discriminant, are situated at infinite values of invariants: 
I4 — > 00 (t — > 0) and 7 6 - * — 00 (* - * 00). 

Hypergeometric integral. Using cq. (22) once again, we obtain 



'3|3 



' S -7/12(l_ s) -H/12 

«s -^r^s „ ^ — + c 2 / 6 



(32/| + 3aJ|)Vi2 



-/ 



ds 
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(32/| + 3^2)7/12 



(36) 



Just as in the case of J 2 |4, it is possible to make the transformation (25), which gives 



4|3 ) (5) = ( J D3| 3 ) 
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We conclude, that cases 2|4 and 3 1 3 literally correspond one to another. From the point of view of discriminant 
theory, this is a remarkable correspondence: discriminant D 2 |4 is a simple two-dimensional discriminant, which 
is well-known and studied, while Z? 3 | 3 is a three-dimensional discriminant, much more complicated and less 
widely known. The study of integral discriminants reveals a parallel between these two cases. 

4 Conclusion 

In this paper we have described the first steps into the study of non-Gaussian averages 
Jn\r(S) = J dx\ . . . dx n e - s ( x i.---» x n) = function of invariants of S 

which are functions only of SL(n) invariants of S, because of inherent SL(n) symmetry. The real motivation 
of this study is generalisation of the well-known Gaussian formula 



Jn\2{S) = j dxi. ..dx n e Si i XiX * = -^=L 



to non-quadratic forms S. If found, such generalisation may immediately have a wide range of applications in 
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statistical physics and quantum field theory. We have worked out several low-dimensional cases and present 
the results in the form of the following table: 



Invariants 



Discriminant D, 



n|r 



Integral discriminant J n i, 



h, h 



It - 64/ 8 



32/| + 37| 



-1/2 



-1/6 



1/4 - i (i/i2) l (5/i2) l fei'i 



Si! 



(1/2), 
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-1/10 _■_ 
4^„ iNI 



-1/2 



-i/4 ~ 1 (l/12) i (5/12) < 
' »! ' (1/2), 



3j|_y 

32/f J 



where 

(a) fe = v = o(o + 1) • • • (a + k - 1) 

i [a) 

As one can see, these series posess a nice structure: they are all hypergeometric, i.c, their coefficients are ratios 
of T-functions. Despite there are too few examples to make far-going conclusions, we conjecture that 

J nV (S) =j dXl ... dXn e-(--) = hypergeometric function of invariants of S 

To support this conjecture, it is necessary at least to calculate a few more examples. As follows from the 
table Fig. 3, the next simplest examples would be J 2 \e, ^3|4 and J 4 | 3 . Also, because of hypergeometric nature 
of integral discriminants, it is interesting to find their q-deformation, as well as the CFT representations of 
these objects, in the spirit of [13]. Especially interesting is the interplay between non-trivial g-deformation 
and action-independence of integral discriminants. Of course, most interesting would be generalization of our 
results to functional integrals, some independent attempts in this direction have already been made [14]. 
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